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We study the electronic structure within a system of phase-decoupled one-dimensional supercon- 
ductors coexisting with stripe spin and charge density wave order. This system has a nodal Fermi 
surface (Fermi arc) in the form of a hole pocket and an antinodal pseudogap. The spectral func- 
tion in the antinodes is approximately particle-hole symmetric contrary to the gapped regions just 
outside the pocket. We find that states at the Fermi energy are extended whereas states near the 
pseudogap energy have localization lengths as short as the inter-stripe spacing. We consider pairing 
which has either local d-wave or s-wave symmetry and find similar results in both cases, consistent 
with the pseudogap being an effect of local pair correlations. We suggest that this state is a stripe 
ordered caricature of the pseudogap phase in underdoped cuprates with coexisting spin-, charge-, 
and pair-density wave correlations. Lastly, we also model a superconducting state which 1) evolves 
smoothly from the pseudogap state, 2) has a signature subgap peak in the density of states, and 3) 
has the coherent pair density concentrated to the nodal region. 

PACS numbers: 74.25.Ha,74.25. Jb,74.72.-h,79.60.-i 



I. INTRODUCTION 

The mysterious normal state pseudogap phase may 
hold the key to understanding the origin of high- 
temperature superconductivity in cuprate superconduc- 
tors. The basic phenomenology of the pseudogap state 
is that of a partial suppression of the density of states 
(DOS) around the Fermi ener gy l eaving only a residue 
of the expected Fermi surface!^] The gap as mapped 
out by angle resolved photoemission (ARPES) has an 
angular dependence which is similar to the d x 2_ y 2 form 
that is realized in the superconducting state with the 
exception of a remaining partial section of Fermi sur- 
face (so-called Fermi arcP) in the "nodal" part of the 
Brillouin zone (BZ)f^ Therefore the pseudogap is "antin- 
odal" in the sense that it has a maximum amplitude cor- 
responding to the maximum amplitude of the supercon- 
ducting gap. Because of these similarites and because of 
the low superfluid density of these materials, it is nat- 
ural to interpret the pseudogap in terms of phase dis- 
ordered superconductivity.^ This scenario agrees with a 
large Nernst signal detected at temperatures above T C P 
However, the fact that the pseudogap extends up to tem- 
peratures well above T c where no evidence of fluctuating 
superconductivity is found seems to rule out this scenario 
as an explanation of the overall pseudogap phenomenon. 
Another natural candidate for the opening of a spectral 
gap is generation of ordered phases, possibly in competi- 
tion with the superconducting state. 

A different interpretation of the pseudogap is that it 
originates from a singlet-triplet gap, or spin gap, due to 
strong nearest neighbor spin correlations related to the 
fact that the system is a doped Mott insulatorP This 
is the resonating valence bond (RVB) scenario.^ Along 
these lines it has also been argued that stripes, which are 
unidirectional spin and charge modulations, may play a 
key role and that the pseudogap would be a signature of 



a singlet pair correlation gap on the hole rich stripesP^ 
Such a correlated state would be similar to the doped 
RVB state but with the distinction of being precondi- 
tioned by the charge inhomogeneity. 

The existence of stripes is well established in 
La 2 _ a; Ba x Cu04 (LBCO) and La 2 - K Sr ;c Cu04 (LSCO) 
[with e.g. Nd co-doping] where both spin and charge 
Bragg peaks associated with the stripe order can be 
clearly identified, presumably enhanced by a structural 
transition in these materials.^ More recently, however, 
evidence is accumulating that also other cuprates exhibit 
similar stripe correlations EEU j n regard the univer- 
sal hour-glass spin excitation spectrum appears partic- 
ularly striking,^ and has a simple explanation within 
stripe modelsPSl Nevertheless, the importance and ubiq- 
uity of stripe correlations remain controversial, partly 
because of the complex disorder present in materials 
like Bi 2 Sr 2 CaCu 2 08+5 (BSCCO) which is, however, very 
amenable to spectroscopic probes like ARPES and scan- 
ning tunneling microscopy/spectroscopy (STM/STS). 

It is well-known that ARPES finds a so-called 
momentum-space dichotomy consisting of coherent (inco- 
herent) nodal (antinodal) quasiparticle states in the un- 
derdoped cupratesP Recent ARPES measurements have 
revealed the presence of two gaps in momentum space: 

1) an antinodal (pseudo)gap which persists above T c , and 

2) a gap in the nodal region which exists only at T < T c 
and exhibits a temperature dependence consistent with 
a standard BCS gapPHU More recently Yang et alW 
found effectively two different types of gaps in the nor- 
mal state at T > T c , a particle-hole symmetric gap in the 
antinodal part of the BZ and a gap without particle-hole 
symmetry in the near nodal region. They also found ev- 
idence for a back bending dispersion in agreement with 
other recent work showing evidence for a nodal pocket in- 
stead of a disconnected arcpSl Consistent with a possible 
pairing origin of the antinodal gap, Kanigel et a/pHiden- 
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tified the existence of particle-hole symmetric Bogoliubov 
bands in the pseudogap state. 

STM experiments have consistently reported signifi- 
cant nano-scale gap modulations in the superconduct- 
ing state with large-gap regions dominant in the un- 
derdoped regime where the pseudogap is known to 
be significant pMUl Using the quasiparticle interference 
(QPI) technique, recent STM measurements on both 
BSCCO and Caa-xNa^CuCbCkj (Na-CCOC) have re- 
vealed an electronic dichotomy similar to that found 
by ARPES: in the low-energ y sect or "conventional" dis- 
persing quantum interferenc e! 28 1 29 1 is detected from elas- 
tic scattering of coherent states near the nodal region, 
whereas incoherent states at energies near the antin- 
odal gap g ive rise to a reduced set of non-dispersive 
peaksPSHSU i n real-space the low-energy local density of 
states (LDOS) is largely homogeneous whereas higher 
bias real-space maps reveal characteristic heterogeneous 
checkerboard or locally stripe-like LDOS modulations. 
This data strongly suggests the presence of extended Bo- 
goliubov states near the Fermi arc an d quasi-localized 
states near the antinodal regions SUED The spatially av- 
eraged DOS evolves smoothly from the supercond ucting 
state into the pseudogap phase above T c I 27 l 38 l 4 U i2l wh ere 
spatial nano-scale (pseudo)gap variations persist pMU Fi- 
nally recent pseudogap low-energy QPI STM studies of 
underdoped BSCCC 38 find a gapless Fermi arc and a 
particle-hole symmetric antinodal gapped region in qual- 
itative agre ement with the above mentioned ARPES 
exp eriment s ! 1 9 1 2 A 

Theoretically, a stripe ordered pseudogap state was 
described in the extreme (and unrealistic) limit of de- 
coupled one-dimensional correlated electron liquids!^ 
Although such an approach can successfully classify 
the possible low temperature phases, it cannot read- 
ily incorporate a more realistic band structure in or- 
der to calculate a spectral function that can be com- 
pared to ARPES and tunneling experiments. However, 
within various approximate schemes, many calculations 
of the spectral distribution of a striped density wave 
have been performedP^ED Recently, these have attracted 
considerable interest in the context of high magnetic 
field transport measurements showing quantum oscilla- 
tions possibly due to severe Fermi surface reconstruc- 
tion from stripe ordering.^^ This interpretation rely 
on electron-like single particle properties of bands which 
have predominant spectral weight in the antinodal region 
which see ms at odds with the existence of an antinodal 
pseudogap! 52 * 53 1 It should also be noted that with realistic 
values of the Coulomb repulsion U, both standard mean- 
field theory and more sophisticated approaches seem to 
indicate that there is no nodal spectral weight but only 
antinodal weight related to states centered on the hole 
rich stripes.^ 

The more unsettling problem with these studies is the 
antinodal spectral weight itself. This is expected to be 
a robust feature, arising from the "stripe bands" which 
are the midgap states induced by the anti-phase domain 



walls in the antiferromagnetic spin texture. For finite 
hole doping the spectral weight from these domain walls 
is antinodal for bond aligned stripes as we consider here, 
at least in mean-field. (For realistic stripe periodicity 
versus hole doping.) Thus, even though there is plenty of 
evidence that stripe correlations may be a generic feature 
of these materials it seems clear that such spin and charge 
density wave correlations will not give rise to an antinodal 
pseudogap. 

In this work we approach the pseudogap problem by 
formulating a non-correlated single-particle caricature of 
a striped pseudogap state by including a BCS pair term 
with finite amplitude only on the hole rich stripes and 
which is decoupled by phase disorder between stripes. 
In Hartree-Fock theory of the Hubbard model a striped 
state is realized by solving self-consistently for a collincar 
spin density wave (SDW) and a charge density wave 
(CDW) with the periodicity of the former being twice 
the latterP It is in this framework that we will intro- 
duce a pairing term acting only in the high hole density 
regions of the CDW, which is meant to mimic the actual 
singlet pair correlations below T* of a strongly correlated 
system. 

Our main result is that such a phcnomcnological model 
with stripe order and phase disordered on-stripe pairing 
is in fact broadly consistent with the salient spectroscopic 
features found by ARPES and STM. There is an ap- 
proximately particle-hole symmetric antinodal spectral 
gap and a nodal region Fermi surface (Fermi arc) in the 
form of a hole pocket. In the near nodal region, just 
outside the pocket, the gap is due to SDW order and 
therefore not particle-hole symmetric. The DOS is sup- 
pressed in accordance with the antinodal gap and within 
an energy scale corresponding to the on-stripe pair po- 
tential. In this pseudogap state the gap is positioned 
at the Fermi level, and the DOS is roughly particle-hole 
symmetric but may have a minimum slightly above the 
Fermi energy due to the band dispersion giving the nodal 
pocket. The low-energy states that make up the Fermi 
surface are extended while states near the pseudogap en- 
ergy are strongly localized on the order of the inter-stripc 
spacing. Finally, we also model a e?-wave superconduct- 
ing state coexisting with the spectral pseudogap and find 
a sub-gap peak in the DOS which is a result of gapping 
the nodal arc, and a coherent pair density that is con- 
centrated to the nodal region. 



II. MODEL 

We consider a tight-binding model in a static field cor- 
responding to unidirectional (striped) SDW order and 
a spatially modulated BCS pair field which is uniform 
along the stripe direction but not, in general, periodic 
in the direction transverse to the stripes. We will con- 
sider only bond centered charge period four stripes (spin 
period eight) and, in order to reduce the number of pa- 
rameters, we do not include an explicit field that couples 
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to the charge density.^ Here we will choose parameters 
such that there is a nodal hole pocket for realistic doping 
levels. In addition there will be antinodal weight in the 
form of open orbits or electron pockets. 

On a square lattice with nearest neighbor hopping t, 
next-nearest neighbor hopping t' , and chemical potential 
jj, we have 



k,a Kj"' 



(1) 



ka 



with £^ = — 2t(cos(k x ) + cos(k y )) — At' cos{k x ) cos(k y ) — /i. 
We will take t = 1 and use t' = -0.3. To generate a SDW 
we include the term 



H S dw = m X! <7(-l) y V(x)n x . 
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x,y,cr 

TO aV 1 C k x ,k y! a C k x -q,k y -x,a- 
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In the following, results are presented for bond-centered 
stripes with V(x mod 8) = (.5, .5, —1, 1, —.5, —.5, 1, —1), 
giving finite V q for q — ±7r/4 and q — ±3n/4 only. The 
variation of the SDW is chosen such that the actual calcu- 
lated spin density comes out with a similar variation. The 
spin potential will generate a CDW which has slightly 
higher hole density on the sites with aligned spin poten- 
tial, which we will refer to as the stripes. 

The main purpose of the present study is to investigate 
the effect on the spectral distribution from local pairing 
that acts only on the stripes. For this purpose we include 
the following BCS term 



Hbcs,cL — Arf 2_j e 
n,{y,y'),j=l,2 



i<t>n J J 
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H.c. = 
H.c, 



where (yy'} is sum over nearest neighbors along the 
stripe, <j> n is an arbitrary phase and Xj <n with j = 1,2 
indicate the left (j — 1) and right (j = 2) legs of the 
bond-centered stripe n, with n = 1, ...,N x /4 and N x 
the system size transverse to the stripe extension. The 
notation j refers to the leg opposite to j on the stripe, 
i.e. j = 1 if j = 2 and j = 2 if j = 1. Note that the pair 
creation and anihilation terms are symmetrized in real 
space coordinates and consequently spin singlet. (Never- 
theless, as follows from symmetry considerations, in the 
presence of the SDW term there will in principle be in- 
duced triplet pair correlations at finite momenta. These 
are a natural consequence of the coexistence of SDW and 
pair correlations!^) I* 1 Eq.pl), 



transform of the local phases for the pairing along x and 
y directions respectively. We will mainly consider the 
problem with quenched uncorrelated random phases for 
which in general all components q of the pairing will be 
non-zero with a periodicity q — > q + 7r/2. The above ex- 
pression for Mq is "d-wave like" in the sense that the same 
term defined uniformly over the system would reduce to 
M q (k x ,k y ) = 6 q fi(cos(k x ) — cos(k y )). Based on the un- 
derlying assumption of pair correlations tied to nearest 
neighbor singlet correlations and strong on-site repulsion 
we expect such a d-wave term to be the most realistic. 
Nevertheless, for comparison we also study "s-wave like" 
local pairing 



H BCS , S =A S £( e ^"4,, n , y 

n,y,j=l,2 



1 ) + H.c. 

H.c. , (4) 



with D y (q) as above. Figure m shows the SDW field to- 
gether with the d-w&ve paired Donds in a stripe unit cell. 

The model H = Hq + Hsdw + Hscs.a with a = d.s 
has long-range superconducting order along the stripe 
direction, which is clearly not realistic even in a stripe 
ordered system. Primarily this assumption is made to be 
able to numerically study large systems in the transverse 
direction but may be partly motivated by the expectation 
of longer range pair correlations along the stripes. It 
constitutes an interesting project to extend the present 
model to include more realistic stripe disorder similar to 
the recent studies of Refs. 

For system size N x the full Hamiltonian with a ran- 
dom phase < <f> n < 2ir on each stripe n = l,N x /4 
corresponds to a 4N X x AN X matrix for each momentum 
k y (including k y + 7r), giving a spectrum of Bogoliubov 
quasiparticles. From this we calculate the (zero temper- 
ature) spectral function A(k,u>) and the corresponding 
DOS p(w) — ^^A(k,uj). We average the spectral func- 
tion over several independent realization of the random 
phases. We present results for magnitude of the SDW 
field to = t/3 and compare the cases without pairing, 
with d-wave pairing = i/4, and with s-wave pairing 
A s = t/3. The magnitude of to is chosen such that for 
relevant doping levels, say x < 20%, the stripe ordered 
system exhibits nodal spectral weight. The values of the 
pairing are chosen to be in qualitative agreement with an 
experimentally measured pseudogap Ao which may be a 
large fraction of t. 



III. PSEUDOGAP STATE 



In the pseudogap phase, the low-energy spectral weight 
is presented in Fig. [2] for the three cases of 1) no pairing 
M q (k x ,k y ) = e iq/2 D x (q) cos(k x + q/2) - D y (q) cos(ky) [Fig. [2^,,b], 2) local d-wave pairing [Fig. [2j;,d], and 3) 

local s-wave pairing [Fig. [2j;,f ] . As seen from Fig. [2^,,b 
the stripe ordered system without any local pair correla- 
tions clearly exhibits antinodal spectral weight as found 



with D x (q) = (1/Ay £^V^-<? 4 ") and D y (q) = 
(V^OEJe^"^ 4 ™' + e l ^"-«( 4n+1 ») is the Fourier 




FIG. 1: (Color online) Unit cell of the SDW with the arrow 
length and direction proportional to the local spin together 
with the local d-wave paired bonds on two stripes with indi- 
vidual (and generally distinct) phases (j> n and 4> n +i 
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FIG. 2: (Color online) Intensity plot of the spectral weight 
at the Fermi energy in the first quadrant of the full BZ (units 
of n I a with a the lattice spacing) for m — i/3 and = — l.Ot 
(left column) or = ~l.lt (right column). Top row (a, b) 
is without pairing, middle row (c, d) with phase disordered 
local d-wave pairing = t/4, and bottom row (e, f) with 
phase disordered local s-wave pairing A s = i/3. Arrows in 
(a) and (c) indicate cuts shown in Fig. [4] A Lorentzian 
broadening with width 77 = O.Olt was used. All plots have 
the same absolute color scheme. Calculations are made for 
a system size N x = 200 with periodic boundary conditions 
and averaged over three realizations of the quenched random 
phases on the charge stripes in (c-f). 
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FIG. 3: (Color online) Spectral weight at the Fermi energy in 
the pseudogap state symmetrized with respect to the stripe 
orientation. The parameters are identical to Fig [2J; and [2ji, 
respectively. The "arc" length grows with doping but is in 
fact restricted by the Bragg planes at k — 7r/4 for the present 
case of period-eight spin stripes. 



previously) 44 1 45 1 However, the local pairing on the stripes 
removes this antinodal weight from the Fermi surface as 
seen from Fig. [2j>f. This is a main result of this paper, 
local on-stripe pairing naturally gives an antinodal gap 
in the spectral function. At the same time the nodal hole 
pocket remains largely unaffected. Results for two differ- 
ent chemical potentials are presented in Fig. [2] [compare 
columns] showing that the nodal hole pocket grows with 
hole doping as expected. The actual hole density for 
these simulations come out to approximately 16% and 
20% respectively with an average hole density variation of 
2% between hole rich (stripes) and hole poor regions (i.e. 
around 10% variation around the mean). The spin den- 
sity, 1(5^ y )\ varies in magnitude between approximately 
0.05 on stripes and 0.13-0.15 between stripes for these 
parameters. 

There is, at the level of the residual Fermi surface, 
no clear distinction between d- and s-wave pairing on 
stripes, indicating, as we will see, that it is the local 
nature of the pairing that is essential. As discussed at 
the end of this section the phase disorder is essential in 
producing the gapless Fermi surface sections, a system 
with ordered phases will in general only have nodal points 
at the Fermi energy. In the rest of this section we analyze 
these results further. 

As mentioned above, more realistic simulations of 
stripe correlations in the cupratc materials need to in- 
corporate the role of quenched disorder and concomitant 
generation of short-ranged stripe domainsP^H One im- 
portant result of such disorder averaging is that the orig- 
inal fourfold rotational symmetry of the crystal lattice is 
restored globally. Within the present approach we can 
obtain a similar though admittedly more rough disorder- 
averaging by x •<-> y symmetrizing our data. Performing 
such a symmetry operation of the spectral weight leads to 
the characteristic k x f-> k y symmetric Fermi arcs shown 
in Fig. [3] As opposed to (71- , it) ordering scenarios such as 
conventional antiferromagnetic order or d-density wave 
order within the present approach the Fermi arc is not 
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FIG. 4: (Color online) Spectral weight along the three cuts 
indicated by arrows in Fig. [JJi mnd^p. The top row (al,a2,a3) 
[bottom row (cl,c2,c3)] corresponds to the parameters used 
in Fig. [2^ [2J;]. The band in (al) [along arrow 1 in Fig. [2^] is 
split by the local pair potential A<j = 0.25t giving a pseudogap 
of similar magnitude in the antinodal region (see also Fig. [7|. 
The cuts in the nodal region (a2,a3,c2,c3) [along arrows 2 
and 3 in Fig. [2|i] are less affected by the pair potential and 
there is clearly (see e.g. c3) a band dispersing through the 
Fermi energy giving rise to the nodal hole pocket. Note the 
broadening as a signature of the pseudogap energy scale in 
near-nodal dispersion (c2). Energy u) in units of t and with 
respect to the Fermi energy. 



symmetric around the AF zone boundary [lines connect- 
ing (±7r, 0) and (0, ±7r)] in agreement w ith re cent ARPES 
measurements in the pseudogap phase DIM] i n principle, 
there is a back side of the hole pocket that might be 
detectable^ 1 but we an ticipate that disorder will further 
reduce its weight .ElMI 

Figure [4] displays several energy versus momentum cuts 
[see arrows in Fig. [2^,,c] of the spectral intensity further 
showing the effect of the disordered pair potential. For 
the antinodal cut in Fig. |4j:l a pairing gap is opened 
with a clear approximate particle-hole symmetry of fi- 
nite spectral weight symmetrically around the Fermi en- 
ergy. This is also evident from fixed momentum cuts (en- 
ergy distribution curve (EDC)) in Fig. Note that 
the minimum in the EDC curves in Fig. [5j;l remains 
at the Fermi level in the antinodal region in agreement 
with ARPES data and contrary to most other ordering 
scenarios for the pseudogap phase with non-zero order- 
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FIG. 5: (Color online) Fixed momentum cuts of spectral 
function intensity from Fig. [4|T-c3. 
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FIG. 6: (Color online) The same nodal dispersion as in Fig. 
[4|i2 and [4]:2 plotted over a larger energy range. One sees a 
small flattening of the low energy dispersion and a broadened 
high energy dispersion in c2 as compared to a2 with an abrupt 
change as a function of energy. 



ing vector and no pairing correlations.^ By contrast, the 
low energy nodal weight is much less affected by the pair- 
ing as seen by comparing e.g. Fig. [4^,3 to Fig. |4j:3; in 
both cases there is a clear signature of a band dispersing 
through the Fermi energy giving rise to the nodal hole 
pocket. This is also seen from the near-nodal and nodal 
EDCs shown in Fig. [5j:2 and[5j:3, respectively. At ener- 
gies below the pseudogap scale there is a signature of the 
disordered pairing in the form of a kink in the dispersion 
and a sudden severe broadening of the spectral function. 
This renormalization of the nodal dispersion in the pseu- 
dogap phase with local pairing is seen more clearly in 
Fig. [6j Interestingly these features of the nodal disper- 
sion are similar to what has been reported by ARPES 
measurements and discussed in the literature mainly in 
terms of bosonic mode-coupling scenarios.^ 

In Fig. [7] we show the spatially averaged DOS in the 
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FIG. 7: Density of states without (dashed) and with (solid) 
local pair potential for the same parameters as in Fig. [2^, and 
[2J:, respectively. The magnitude of the local pair potential 
Ad = 0.25t is clearly reflected in suppression of the low-energy 
DOS below a pseudogap scale Ao « 0.25t and is also seen to 
correspond to the splitting of the antinodal band (Fig. |4j;f ) 
around the Fermi energy. 
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FIG. 8: (Color online) Nodal (solid), antinodal (dashed), 
and total (dotted) low-energy DOS. We find an approximately 
particle-hole symmetric antinodal pseudogap and a DOS con- 
sistent with a band dispersing through the Fermi energy in 
the nodal region. The nodal weight has a minimum at the 
edge of the nodal band where the pocket closes. Inset shows 
the spectral weight (c.f. Fig. [5J;) at the Fermi energy and the 
shaded (white) regions of the BZ are integrated over to get 
the antinodal (nodal) DOS. 

case with (solid) and without (dashed) pairing correla- 
tions on the stripes. One can identify the pseudogap 
scale A ~ by the suppressed spectral weight be- 
low this energy. The magnitude of the pseudogap in the 
DOS clearly corresponds to the main antinodal gap in 
the spectral function. Even though the overall shape of 
the calculated spatially averaged DOS in the pseudogap 
phase consists of a suppression of spectral weight at the 
Fermi level, an interesting feature is that the minimum 
is not necessarily tied to the Fermi energy but shifted 
slightly to positive bias u> > 0. This is a manifestation 
of the fact that the nodal band disperses through the 
Fermi energy and consequently will not have a minimum 
at ui = 0. This, we suggest, is a telltale sign of the quasi- 
particle nature and band dispersion of the nodal states. 
Interestingly, recent STM experiments may already have 
meas ured such a shift to positive bias in the pseudogap 
statePCH In our simulations we can study nodal and 
antinodal contributions to the DOS directly by integrat- 
ing the spectral weight over only nodal or antinodal re- 
gions. For this purpose we define 

7T 7T — fco 

4 E E A ^)i (5) 

Pantinodal{^) = p{u) ~ Pnodal(u) , (6) 

where kg is some appropriate division between nodal and 
antinodal parts of the BZ. Such a division of the low- 
energy spectral weight is natural based on the discon- 
nected nature of the Fermi surface in the stripe ordered 
system as seen e.g. in Fig. [2^, and [2]:. Figure [8] displays 
this division of the DOS for the same system with phase 
disordered d-wave pairing as shown by the solid line in 




FIG. 9: (Color online) One dimensional (direction transverse 
to stripes) IPR as defined in Eq. Q versus energy of quasi- 
particle excitations E-y for the same parameters as in Fig. [2J; 
with system size N x — 200 (small points) N x — 320 (large 
points). A large fraction of the states near the pseudogap 
energy Ao ~ Ad = 0.25i are strongly localized in contrast to 
states near the Fermi energy. Inset is the corresponding IPR 
for the same parameters but with the same pairing phase on 
all stripes resulting in all states being extended Bloch states. 

Fig. [7] (using k — 0.37r). The nodal DOS is consistent 
with a band dispersing through the Fermi energy with 
the minimum (uj ss 0.05) corresponding to closing of the 
hole pocket (see Fig. Ilb3). By contrast, the antinodal 
DOS has the property of an approximately particle-hole 
symmetric pairing gap. From this division it is evident 
why the total DOS may have a minimum shifted to posi- 
tive energies. As a point of principle this is an important 
example, however in general we do find that the DOS 
in the pseudogap phase for other parameter values may 
have a minimum at u> = 0. Even though there is a band 
dispersing through the Fermi energy it may have a sup- 
pressed DOS because it is affected by the pair potential. 

To further understand the nature of the low-energy 
states within this pseudogapped model system we study 
to what extent the quasiparticle excitations are localized 
as an effect of the disordered pair potential. To this end 
we calculate a one dimensional inverse participation ratio 
(IPR) since we can only have localization transverse to 
the stripe extension in the present model. Because of the 
local superconducting order, the quasiparticle excitations 
are Bogoliubons which give a contribution to the LDOS 
at both positive and negative energies. To get a measure 
of the extent of a state in the transverse (x) direction we 
add up these contributions and average along the stripe 
direction (y): 

IPR 7 = C£lJ2(^y^ 2 + ((cL,^ 1 )) 2 ] 2 )- 1 • (7) 

x y,a 

Here the average (•) is taken with respect to the ground 
state and 7^ is the creation operator of a quasiparticle ex- 
citation with energy E 1 . A completely delocalized quasi- 
particle with equal weight on all sites (actually chains 
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along y) will have IPR = N. x , with N x the system size, 
whereas a quasiparticle which only has spectral weight 
on a single site will have IPR = 1. Figure [9] is a plot of 
IPR as a function of quasiparticle energy of all states in 
an energy window uj < 0.5t for system sizes N x = 200 
and N x — 320. As seen, states around the pseudogap 
energy = 0.25£ are most localized, some with an IPR 
as small as the stripe spacing and independent of sys- 
tem size, indicating that the states are basically local- 
ized on a single stripe (we expect localization length £ ~ 
IPR/2). Strikingly, the very lowest energy states (which 
make up the nodal hole pocket) all have a much larger 
IPR (here > 70) , and are clearly only weakly affected by 
the stripes and the disordered pair potential, in fact most 
of these states scale with system size. This is fully con- 
sistent with the band-like electronic quasiparticle nature 
of these states as discussed above. 

The transverse localization of the antinodal states at 
the pseudogap energy also makes it clear why the on- 
stripe pairing acts so strongly on these states; these states 
have most of their spectral weight on the stripes. This 
property of the antinodal stripe states was emphasized 
earlier for the case of diagonal disorder and related to 
the narrow band width in the transverse stripe direction 
and c orresp onding quasi-one dimensional nature of these 
statesP^^ID xhe role of the phase disorder in the present 
calculations is thus twofold, 1) it acts to localize antin- 
odal states on stripes thus making these more susceptible 
to the local pairing, and 2) it causes the more extended 
nodal states which have most of the spectral weight be- 
tween stripes to see only a very small spatially averaged 
pair potential, thus making these states relatively insen- 
sitive to the local pairing. 



IV. SUPERCONDUCTING STATE 

So far we have considered only phase disordered pair- 
ing, suggesting that this is a caricature of a correlated 
pseudogap state with local pair correlations in the form 
of an on-stripe spin gap. The most obvious extension 
of this model to a state with long-range superconduct- 
ing order would be to lock the phases between stripes. 
Such a state would correspond to the mean-field theory 
of a stripe ordered superconductor with a periodically 
modulated superconducting order parameter, with finite 
components at q = and q = ±7r/2 (not q — tt for bond- 
centered stripes). We will not study this state here and 
only remark that in contrast to the phase disordered pair- 
ing, this state (for the d-w&ve case) has a Fermi surface 
consisting of points at the nodes of the gap function. For 
the particular case considered here (with pairing only 
on the bond centered stripes) we find a nodal line at 
2cos(fc y ) — cos(k x ) — thus resulting in a point node 
that is shifted slightly away from that of an actual d-wave 
superconductor with intact fourfold rotational symmetry. 

Another alternative, similar to that suggested in 
the context of stripe ordered x = 1/8 doped 




FIG. 10: (Color online) Density of states (curves off-set for 
clarity) for different magnitudes of homogeneous d-wave pair 
potentials A h = (solid), A h = 0.05t (dashed), A h = O.lt 
(dotted). The pseudogap energy is Ao and the energy Ai 
signifies the peak from the gapped nodal pocket in the super- 
conducting state. The other parameters are the same as in 
Fig. fifr: m = t/3, A d = 0.25t, and fj, = -t. The thin line 
shows as comparison the DOS for a system with only stripe 
order (m = t/3, Ad = 0, fi = —t) and homogeneous d-wave 
potential A^ = 0.15t. 



Lai. 6 _ x Ndo.4Sr !C Cu04 La 2 ^ ;E Ba a; Cu04, lsZMI would be to 
lock the phases with a 7r phase shift between neighbor- 
ing stripes, resulting in a striped superconductor without 
a q = pairing component, and only finite momentum 
pairing. Because there is no zero momentum pairing such 
models have an extended Fermi surface similar to what 
we find here for the phase disordered pairingPsl 

In the following, we suggest a different model for the 
superconducting state which, in addition to the phase 
disordered pairing on stripes, also includes a uniform d- 
wave pair potential. The motivation for studying such a 
state comes from considering a glassy stripe state with 
static but only short range stripe order. In the supercon- 
ducting state we expect there may be finite momentum 
pair correlations developing due to the local stripe order, 
but since phase coherence is established between stripe 
patches with different orientations only the q = compo- 
nent may order. Thus in addition to the original model 
with phase disordered on-stripe pairing we include ex- 
plicitly the homogeneous component 



H B cs,h = A h ^[cos(fc y ) 
fe 



cos(M]4 t c U + H - c - ( g ) 



This homogeneous pairing opens a <i-wave gap on the 
Fermi arc. This is the actual gap of the superconduct- 
ing condensate which will disappear in a standard mean- 
field fashion as temperature is increased.^ As a signa- 
ture of the gapped pocket a sub-gap peak appears in 
the DOS as seen from Fig. 10 The origin and po- 



sition of this peak is set by the transition of the con- 
tours of c onsta nt energy from separated "bananas" at 
low energjlM22l ^ a pocket at higher energy. This is sim- 



ilar to what has been discussed recently in the case of 
coexisting d-wave superconductivity and commensurate 
antiferromagnetismP^ The energy of the subgap peak, 
Ai, depends on band-structure and is not simply linearly 
related to the magnitude of Aft as seen in Fig. [l0| which 
shows the evolution of the DOS with increasing magni- 
tude Aft,. The pseudogap energy scale, A , that follows 
from the phase disordered pairing is maintained even in 
the presence of the homogeneous term as long as the lat- 
ter is smaller in magnitude. In contrast, a stripe ordered 
state with only a homogeneous pair potential has strong 
coherence peaks at twice the gap magnitude in addition 
to the smaller peak from the gapped nodal pocket as seen 
from the thin line in Fig. [To} 

Figure [IT] shows the DOS spatially averaged over sites 
on (off) the charge stripes only, i.e the solid curve displays 
the DOS spatially averaged over sites that contain pseu- 
dogap bonds (see Fig. [TJ. One sees a low-energy "uni- 
versal" LDOS consistent with homogeneous real-space 
STS spectra at low (sub-gap) bias, but deviations in the 
LDOS at higher energies resulting in spatially modulated 
(in this case stripy) LDOS real-space maps. The Fourier 
transform of the LDOS at a q vector corresponding to the 
charge order, q = 2tt/4, is proportional to the difference 
between the curves shown in Fig. and hence roughly 
particle-hole symmetric, even in the pseudogap state. 
Interestingly, a plot of the actual zero momentum pair 
shown in Fig. 
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reveals that this 



quantity is dominated by weight in the nodal region on 
the hole pocket even though we consider a d-wave gap 
with larger magnitude in the antinodal regions. The rea- 
son for this is of course the coexisting antinodal pseudo- 
gap. This is consistent with a recent analysis of ARPES 
data by Kondo et aZ.^in which it was found that most of 
the coherent weight in the superconducting state is near 
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FIG. 12: (Color online) Zero mometum component of the 
pair den sity for parameters identical to the dashed curve in 
Fig. [To] m = i/3, A h = 0.05t A d = 0.25i, and ^ = -t. 
The shading indicates the deviation from zero. The weight is 
along the tight-binding Fermi surface with highest intensity 
on the nodal pocket of the pseudogap state. The real part 
is plotted, there is also a smaller imaginary part (not shown) 
from the phase disordered on-stripe pairing. 



the node. Similar conclusions about a coherent smaller 
gap below T c on top of the pseudogap have al so been 
reached in earlier STM, ARPES and ^SR work! 16 ' 27 ' 75 ! 
Such a two gap phenomenology is naturally suggestive of 
a competition between a pseudogap state of non-pairing 
origin and a more traditional superconducting state. In 
the present model where the pseudogap is a pairing gap 
we have in mind an alternative scenario in which the pair- 
ing arises on the stripes which then infects the extended 
nodal states where the actual phase coherent pair den- 
sity is formed. Although this mechanism is not contained 
in the model we may speculate that the antinodal spec- 
tral weight cannot couple coherently, possibly because 
the pseudogapped states are too localized or because the 
corresponding local pair density has most of the weight 
in finite momentum components. 



V. DISCUSSION 



FIG. 11: Local density of states in the superconducting state 
on (solid line) and off (dashed line) a charge stripe. The 
parameters used for this curve are similar to the dotted line 
in Fig. [To] The overall collapse of these two curves at low 
energy results in real-space homogeneity at low tunneling bias 
whereas the difference in LDOS at high energy will show up 
as stripes in the real-space STM field-of-view. 



The model proposed in this paper does not fall into the 
conventional separation of scenarios for the pseudogap 
state: either preformed pairs or competing order. Our 
scenario contains both an order (SDW and CDW stripes) 
and concomitant local singlet pairing states which give 
rise to the spectroscopic antinodal pseudogap. Surpris- 
ingly, however, the Fermi arc on the nodal hole pocket 
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is unaffected by this pseudo-gap mechanism irrespective 
of whether the local pairing symmetry is of s- or d- 
wave type. For the ordered SDW stripes studies here, 
the Fermi arc is the front side of the nodal hole pocket. 
However, in a more realistic short-range glassy p hase t he 
pocket is wiped out and replaced by a single arc! 51 * 64 ! 

As mentioned above, the phase disorder has negligible 
effects on the extended nodal states since the effective 
pairing potential averages to zero for these states. One 
can get additional understanding of the role of phase dis- 
order by comparing to a stripe state with non-zero pair- 
ing and fixed phase [i.e. (p n — for all n in Eq. (3)]. 
In such a phase-ordered stripe state, the pairing leads to 
nodal points at the Fermi energy, and from this respect 
the phase disorder is crucial for generating a Fermi arc. 
The origin of this increased spectral weight at the Fermi 
level can be traced to the fact that significant phase dif- 
ferences (e.g. 7r-phase shifts) between d-wave supercon- 
ducting regions separated by an tiferromagnetic regions 
generates low-energy states ESE3 

A recent study of thermal phase fluctuations of a d- 
wave superconductor also found the presence of a Fermi 
arc.^ In that work too, the phase disorder is crucial for 
generating the arc, but contrary to the present approach, 
no pocket exists in the nodal region from SDW ordering. 
Even though the pocket could be destroyed by disorder 
one may be able to distinguish these scenarios for the 
Fermi arc by searching for particle-hole symmetry in the 
near nodal region, a property that should be absent in 
that region of momentum space within the approach pre- 
sented here. (The latter in agreement with the recent 
analysis of ARPES measurements by Yang et a/P^) 

Experimentally it has also been found that the Fermi 
arc length scales with T/T*^ Even though extension 
of the present model to a self consistent study of finite 
T effects is beyond the scope of this paper, we remark 
that a model of granular antiferromagnetic and d-wave 
superconducting islands found that the phase disorder 
can indeed generate a T-dependence of the Fermi arcpl 
similarly to the work by Altman et al^ For the present 
model it is natural for the arc length to vary inversely 
to the magnitude of the SDW potential which in mean 
field is expected to grow with decreasing temperature. 
Whether or not such a variation can quantitatively re- 
produce the experimentally observed variations of arc 
length with temperature remains to be explored in more 
detailP 

Because the pseudogap proposed here remains at the 
Fermi level (modulo a possible small shift) it explains the 
smooth evolution of the spatially averaged DOS as a func- 
tion of temperature from the s uperconducting state into 
the pseudogap stateP^EHEIEU As we have shown, nodal 
states are extended whereas antinodal states are localized 
effectively on individual stripes in overall agreement with 
recent conclusions from e.g. tunneling spectroscopy.^ 
The nodal band can cause a slight shift of the DOS min- 
imum to positive bias in the pseudogap state similar to 
what has been seen by recent STM measurements at 



T > T c .1221221 From the ARPES evidence of a quasipar- 
ticle band dispersion through the Fermi level together 
with an antinodal gap, this shift is something we expect 
on general grounds, and should be explored further by 
future tunneling experiments. 

The superconducting q — condensate at T < T c 
essentially "lives on top of the pseudogap phase and 
its main spectroscopic effect is to gap the states at 
the Fermi arc. The superconducting gap gives rise to 
a sub-gap structure at T < T c within the pseudogap 
DOS. Interestingly, sub-gap features in the spatially av- 
eraged DOS which disappear above T c have been recently 
pointed out by s everal STM experiments of underdoped 
cuprates! 27 * 81 * 82 * A very interesting aspect of the present 
model for the superconducting state is that it naturally 
incorporates the existence of spatially varying local pair- 
ing amplitudes. This aspect was previously proven suc- 
cessful in describing salient features of th e LDO S gap 
modulations measured by STM at T < T c \ P 1221 Specif- 
ically, it reproduces the presen ce of sharp coherence 
peaks in small gap regions-^I^U and the g ranul ar tran- 
sition through T c observed by Gomes et ai! 26 * 73 * It is an 
interesting future study to extend the present model to 
include more realistic disorder configurations and deter- 
mine whether the spatially resolved LDOS is in further 
agreement with STM data. In the superconducting state 
we have seen that the low-energy LDOS is roughly con- 
stant (in space) within this approach but it is known that 
correlations which penalize charge fluc tuatio ns will fur- 
ther stabilize nodal LDOS universality! 83 * 84 * Also it will 
be interesting to see whether the present approach can 
reproduce the fc-space dichotomy seen by Fourier trans- 
formed STM maps between the low- and high-energy 
momenta.^ In addition an essential aspect of the model 
is that the antinodal "stripe" states are localized by dis- 
order in the phase of the local pair potential. An interest- 
ing question is whether similar localization may perhaps 
more realistically be caused by stripe density wave disor- 
der. 

Lastly, since the present model assumes SDW order it 
includes per construction local moments which, as shown 
recently, are present even in the magnetic response of 
the BSCCO materials.^ The overall hour-glass neutron 
response is therefore fully compatible with the present 
model even though one may have to include glassy dis- 
order and/or soft fluctuations to reproduce details of the 
measured low-energy magnetic fluctuations. 



VI. CONCLUSIONS 

We presented a phenomenological model for the pseu- 
dogap state consisting of stripe spin- and charge density 
wave order with phase decoupled on-stripe singlet pair- 
ing. This phase is characterized by extended states giving 
rise to a nodal hole pocket and gapped antinodal states 
that are localized transverse to the stripes. The spectral 
function exhibits particle-hole symmetry characteristic of 



10 



a pairing gap in the antinodal region whereas in the near 
nodal region just outside the hole pocket the gap is due 
to stripe order and does not display particle-hole symme- 
try. The nodal dispersion has the appearance of a normal 
metal at the Fermi level but contains a kink and broad- 
ening below the pseudogap energy scale. The spatially 
averaged DOS has a gap at the Fermi level and is also 
roughly particle-hole symmetric in the pseudogap state 
even though a small shift arising from the nodal band 
may be observable. Finally we also discussed a super- 
conducting state coexisting with the pseudogap where 
the nodal pocket becomes gapped. Here we find a char- 
acteristic sub-gap peak in the DOS and a coherent pair 
density with most of the weight in the nodal region. 
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Note added. After the submission of this work a new 
ARPES study of the pseudogap state in underdoped 
Bi2201 appeared.^ An extended temperature regime 
ranging from above the pseudogap transition T* to below 
the superconducting transition was studied and found a 
conspicuous shift of the antinodal gap minimum momen- 
tum in the pseudogap state compared to kp above T* as 
well as an unexpected additional spectral broadening in 
the superconducting state. The observed shift suggests 
that the antinodal pseudogap is not (exclusively) a pair- 
ing gap. In light of this new ARPES data we realized that 
these qualitiative features were already present in the 
model considered in this paper. We thus present an ad- 
ditional Fig. ( 13 1 using our earlier calculations showing 
a comparison between antinodal spectral weigth cuts for 
the bare band structure (modeling the region T > T*), 
the pseudogap state, and superconducting state. There 
is a shift of the putative (turning point of coherent 
dispersion) in the pseudogap state compared to kp from 
the bare band structure, which is due to the stripe or- 
der. There is also a broadening of the spectral weight 
in the superconducting state, which is related to the in- 
terplay between phase coherent zero momentum pairing 
and the local incoherent on-stripe pairing. The details 
of this broadening and the generics of the stripe induced 
shift of kp will be studied in more detail in subsequent 
work. 
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